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Abstract 

The temperature phase transition in scalar <f> 4 (x) field theory with spon- 
taneous symmetry breaking is investigated in a partly resummed pertur- 
bative approach. The second Legendre transform is used and the resulting 
gap equation is considered in the extrema of the free energy functional. 
It is found that the phase transition is of first order in the super daisy as 
well as in a certain beyond super daisy resummations. There no unwanted 
imaginary parts in the free energy are found but a loss of the smallness of 
the effective expansion parameter near the phase transition temperature is 
found in both cases. This means an insufficiency of the resummations or 
a deficit of the perturbative approach. 



1 Introduction 

The question on the order of the phase transition when a spontaneously broken 
symmetry is restored by raising temperature is a topic of considerable interest. 
The present status of the problem in scalar 4 -theory with spontaneous symmetry 
breaking due to the mass term with 'wrong sign' is characterized by the general 
opinion that it is a second order transition, see e.g. the standard textbooks [[I], 
[2|. m . This opinion is based on non-perturbative methods like lattice calculations 
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and the average action resp. flow equation method [HI [J, however still lacks a 
formal, final proof. 

In the perturbative approach numerous authors found a first order phase 
transition. Here the problem is in the infrared divergencies appearing there. 
A resummation of the relevant diagrams is necessary. This was known already 
since the pioneering papers in the field, e.g., ||. The general outcome is that the 
perturbative methods do not give reliable results near the phase transition for 
two reasons: 

1. Unwanted imaginary parts, for instance in the nontrivial minimum of the 
free energy (resp. the effective potential) appear indicating a instability 
or insufficient resummation. So for example, in the paper |7j a first order 
transition had been found on the level of super daisy graphs resummed. 
But then the authors included higher loop graphs and found an unwanted 
imaginary part. 

2. Near the phase transition, the effective expansion parameter of the re- 
summed perturbative series may become of order one making the contri- 
butions of higher order graphs unpredictable. For instance, it cannot be 
excluded that their inclusion may change the order of the phase transition. 

In the present paper we consider the scalar 4 -theory in (3+1) dimensions in 
the perturbative approach. We perform a resummation of the perturbative series 
using the second Legendre transform and go beyond the super daisy resummation. 
We refine our previous results || obtained by a simpler resummation method 
which is restricted to the level of including the super daisy diagrams only. The 
disappearance of all unwanted imaginary parts was shown there and, thus, is 
repeated in Sec. [| We note that this result implies that the imaginary part 
found in |7[] from higher loop graphs is questionable. 

Below, in Sec. [|, we calculate a number of contributions beyond the super 
daisy level. We show that, near the phase transition, they are, as expected, of the 
same order as the leading contributions. They turn out not to change the order 
of the phase transition but to affect its quantitative characteristics only making 
the transition a bit stronger first order. 

To explain the controversy results on the order of the phase transition one 
could argue that the higher loop graphs not included in the perturbative approach 
might sum up to contributions that weaken the first order character until it turns 
into a second order one when all graphs are summed up finally. Our result, 
however, points in the opposite direction. 

As technical tools we use the imaginary time formalism and functional meth- 
ods. The most important tool is the use of the second Legendre transform^ 
together with the necessary condition for an extremum of the free energy. 

1 Is known since the beginning of the sixties (j^, [l(J and can be found in textbooks [jllj, see 
also |2)). 
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We consider the scalar 4 -theory in (3+1) dimensions with the 'wrong sign' in 
the mass term which takes after symmetry breaking by shifting the field <j>(x) — > 
v + 4>(x) the form 

s[<t>] = ^ 2 --/ (1) 

+ J dx ^0(x)K^0(x) - Xv(f)(x) 3 - ^0(x) 4 + vim 2 - \v 2 )<\>{x)^ 

with K M = □ — fi 2 and the tree mass p 2 = —m 2 + 3Af 2 . The line A (free 
propagator) in the corresponding graphs is the inverse of K: A = — K -1 . The 
quantity to be calculated is the free energy F. By means of W = —TF it is 
related to the vacuum Green functions 



Z = j D(j)e b , (2) 

where T is the temperature and its connected part reads 

W = \ogZ. (3) 
After performing the second Legendre transform the representation 

W = S[0}+ l -Tr log - ^TrA^p + W 2 [p] (4) 

emerges, where all graphs are to be taken with the line P (instead of A) which is 
subject to the Schwinger-Dyson (SD) equation 

p-\ P ) = A- i -m(p), (5) 

where 



are the self energy graphs with no propagator insertions, see, e.g., [0 for details, 
is the sum of all two particle irreducible (2PI) graphs . 
In the imaginary time formalism the operation 'Tr' is given by 

with the momentum is p = (2nTl,p). We indicate the dependence on a functional 
argument P by square brackets and on the momentum p by round ones (P(p) is 
the corresponding Fourier transform), 5 is the variational derivative. 

Turning to perturbation theory it is useful to start from W%. The first contri- 
bution to its perturbative expansion is the graph looking like a 'eight', 

r(l) 



wr = on- (8) 
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By means of the formulas given above it generates the free energy with all super 
daisy diagrams summed up. In this way the results of our previous paper || are 
reproduced. In the present paper we go beyond the super daisy level by including 
the following (higher loop) graphs 

^ = 1 8 + h ® + s : 2> + £ f a + E a • (9) 

n>3 n>3 

where the dots symbolize further subgraphs to be inserted so that the number of 
vertices is n. This selection of graphs is motivated by the 1/N expansion where 
it includes all leading and first non-leading contributions. 



2 The SD equation and the extrema of the free 
energy 

The problem to be considered with respect to the phase transition is whether 
the free energy exhibits a nontrivial (i.e., at v ^ 0) maximum and minimum for 
some v > at finite T. The common approach is to solve the SD equation (|5|) 
in some approximation and to insert its solution into the free energy (§) with 
a subsequent investigation of its dependence on the condensate v. Instead, we 
suggest to consider the necessary condition for the extremum of the free energy 

dW 

a^ = ° < 10 > 

together with the Schwinger-Dyson equation (|5]). In this way one variable can be 
eliminated and we obtain a simpler expression. In fact there are more simplifica- 
tions behind. Away from the extrema of the free energy there are tadpole graphs 
generated by the linear term in the action ([I]) which must be taken into account. 
They disappear once Eq. ( |TUD holds. 

We obtain from Eq. fllOD using (£|) and A -1 = — K^ 1 = p 2 — m 2 + 3Xv 2 the 
condition 

dW m 2 A 2 3 dp (1 , 1 SW 2 \ dW 2 

which by means of fl5|) reduces^ to 

\ v 2 = m 2 - 3AS + 2— -f 12 

ov 2 



with the notation 



] = Tr/3. (13) 



2 This is related to the known fact that the SD equation (^) is a stationary point of W 
considered as a functional of /3, see e.g. p3. 
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A representation of So is given in the Appendix. 

Now the SD equation fl5|) can be rewritten eliminating v using (|T2|) to become 

. dW 2 
dv 2 



/T 1 = p 2 + 2m 2 - 9A£ + 6^ - E[/3](p) . (14) 



It holds in the extrema of the free energy and allows directly to determine number 
and kind of the extrema. 

Being interested in the infrared behavior at high temperature we make the 
ansatz 

(3~\p) =p 2 + M 2 (15) 
and obtain from Eq. flT3|) the gap equation in the extremum 

M 2 = 2m 2 - 9A£ + 6^ - E[/3](0) , (16) 

where the ansatz (Jig) is assumed to be used for (3 in the right hand side. 

To proceed we need now to make the approximation defined by Eq. (|^) on 
the level of the 2PI functional. The analytic expression for W 2 corresponding 
to Eq. (D) can be obtained from the graphical representation by inserting the 
line P and the vertex factors — 6Xv resp. — 6A for the three resp. four vertex. It 
reads 

W£ h) = -^AS 2 + 3A Vr + \ 2 D , (17) 
where we introduced the notations 



and 



as well as 



(19) 



Z 1 (p)=Tr q P(q)f3(p + q)= -<> . (20) 

Note that the infinite set of diagrams in Eq. @ can be reproduced by a formal 
expansion of (|T7|) with respect to powers of A. 

From these formulas we obtain by means of (||) the self energy in the form 

E[ /3 ](p) = -3A So + 6 AViL + 2A^, (21) 

where we used obvious relations like 8f3 (p) / 5(3{q) = S(p—q) and 5T*i [(3] (p) / 5f3(q) = 

2/3{p + q). With gja = 3A 2 r and Eq. (||) we obtain the SD equation in the 
extremum 



/T 1 = p 2 + 2 (m 2 - 3AS + 6A 2 r) ^1 - 3A ^y) + 6A2r - 2 ^T^ . ( 22 ) 



8D 
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and from fllB] ) the gap equation 

M 2 = 2f (m 2 - 3A£ + 6A 2 r) (23) 
in the given approximation. Here we used the relation 

6A2r - 2X m - ^ - «> + *» E ^lrli|y (24) 

showing that this contribution to Eq. (^) vanishes for p = 0. The quantity T is 
given by 

5(3(0) l + 3AEi(0) q K J (l + 3A£i(g)) 2 ' v ; 



where T fll8|) was used. 

Note that Eq. (^) does not contain the condensate v so that it is an equation 
to determine the mass M in the extrema of the free energy as a function of A, T 
and m. 



3 Symmetry restoration in super daisy approx- 
imation 

Now, in order to establish the connection to our previous results || and to make 
the subsequent calculations more transparent we restrict ourselves for a moment 
to the level of super daisy graphs summed up only. This is obtained by setting 
r and D in the formula of the preceding section equal to zero. So, with 



W^> = g Q = (26) 

we obtain the gap equation in the extremum in the form 

M 2 = 2m 2 - 6AS (27) 

(cf. Eq. (37) in 0). For the propagator = p 2 + M 2 is exact because S does 
not depend on the momentum p. Equation (|l|) for the condensate reduces to 

\ v 2 = m 2 - 3AS (28) 

so that we can note 



, 9 M 2 

Xv 2 = — (29) 



in the given approximation. 
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The solution of equation ( p?D can be investigated graphically by drawing 
both sides of the equation as function of M in one plot. A detailed analysis had 
been given in |J showing clearly the behavior expected for a first order phase 
transition. For instance, in a certain temperature region the existence of two 
solutions corresponding to a maximum and a minimum of the free energy was 
found. For small coupling A, which corresponds to high temperatures of the phase 
transition T ~ A -1 / 2 , the function S in the gap equation can be expanded in 
powers of T/M, see Eq. (|55|) in the appendix for details. Keeping the relevant 
orders the gap equation turns into 

„ „ (T 2 MT m 2 M 2 / (4vrT) 2 \\ , x 
M 2 = 2m 2 -6A(-- — + — + — (log^-^-2 7 )) (30) 



and has the solutions 

^ Min ~ 

Max 

where the notations 
and 




f3XT\ 1 „ . XT 2 
I » + 2™%- — , (31) 

'3 2 



5 3 \ { 4tt ) 5 3 



<h - 1-14 ( 32 ) 



are used. Here, M Max resp. M Min (and v by means of Eq. (p9|) ) correspond to 

the minimum resp. to the maximum of the free energy at v ^ 0. The lower 

spinodal temperature T_, i.e., the temperature at which the maximum appears 

at v = 0, follows from Muax = to be T_ = ^v^i an d the upper spinodal 

temperature T+, i.e., the temperature where the minimum and the maximum 

merge and disappear for T > T + follows from (|3l|) to be T + = ^= , L_ 5 \x ■ For 

Y 87^63 

small A, in the first nontrivial order, they are 

2m / 3A \ 2m f 3A \ 



^ = 1 + 7^5 (35) 



and their ratio 

T_ 167T 2 

is larger than one. We note that this does not depend on the logarithm of T 
coming in with 63. From Eq. ( |3lD it is easy to calculate the masses in T_ and in 

3a/A 

M„. , = m, M„ 1 =0, 

o /T 

M M . , = M M , = m. (36) 
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The temperature T c of the phase transition where the free energy is the same in 
the trivial (i.e., at v — 0) and in the nontrivial minimum can be calculated in the 
given approximation quite easy. For this one needs the gap equation (|5|) at v = 0, 
/T 1 = p 2 - m 2 + 3AE . Note, that the self energy in the given approximation 

^(P) = 2 tJ-T f-7 AS o) = -3AS (37) 



5(3{p) V 4 

is independent of p and the ansatz /3 _1 (p) = p 2 + Mq is exact again, where M is 
the mass as solution of the gap equation at v = 0. In this way the gap equation 
there reads Mq = —m 2 + 3A£o- Taking the approximation fl55|) for So we obtain 

in a way similar to Eq. © its solution: M = =ff + ^(f£g) 2 - - 

Now T c can be obtained from equating the free energy at v = and in the 
minimum: 

m m f Wo = W(u^0,M^M o ) 

<U Wf > with = ^^-^^^ (38 ) 

where are the connected Green functions (f|). We rewrite them using A -1 /? = 
1 — /3T, following from the Schwinger-Dyson equation (ffy and drop the infinite 
constant Trl appearing thereby. Again, in the given approximation, (BBT) , using 



(37) we obtain 



w m = _™L V * + V _ 3 AS 2 + I Trln/3 . (39) 

2 4 4 2 v ; 

Now, Wq resp. Wj 1 ^ are obtained by inserting the corresponding values of v 
and M. After that, Eq. ( j38[) can be approximated in the lowest nontrivial order 
in A (thereby, in |Trln/3 all orders of A shown in (|53| ) have to be taken). The 
temperature T c is in between T_ and T + , closer to T + . With T c = T + (l — Ar) we 
obtain r = 0.0028 as numerical solution of the emerging higher order algebraic 
equation. 

So, on this stage of resummation, i.e., after summation of all super daisy 
diagrams which is generated by taking for W2 the the representation ( p6f) , 
we see clearly a first order phase transition. In the nontrivial minimum the mass 
Mmiu fl3lD is positive and of order y/\: 

3 3 

— y/Xm = M Mi n\ T < M Min < -V\m = M Min \ . (40) 

This mass has to be inserted into all higher loop graphs which are not included 
in the given approximation. As the theory is Euclidean and the propagator is 
j3(p) = l/(p 2 + M 2 ) a consequence is the absence of unwanted imaginary parts. 
In this way the first problem mentioned in the introduction is settled. 
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Consider now the second. A simple counting of the expansion parameter 
of the perturbative series after all super daisy diagrams are summed up goes 
as follows. Let V, C resp. L be the number of vertices, loops resp. line of 
a graph. The vertex factors are proportional to A (0 4 -vertices) resp. Xv (0 3 - 
vertices). By means of formulas and fl36|) the condensate is of order one, 
v ~ m, so that all vertex factors are of order A. The summation/integration 
measure associated each loop is TY^iZ-oof dp/(27r) 3 (|7|) and the factors of the 
lines are l/((2vr/T) 2 + ft + M 2 ). For high T (~ A~ 1/2 ) (||) the graphs can be 
approximated by taking the zeroth Matsubara frequency (I = 0). Than, doing 
a rescaling p — > pM, a factor A = \ V T C M 3C ~ 2L appears in front of each graph 
which itself is then a number of order one. By means of C = L — V + 1 and 
M ~ vA we obtain A ~ A so that it ceases to depend on the order of the 
perturbative expansion. This is a loss of the smallness of the effective expansion 
parameter. Although we know that all graphs are finite and of order one, we do 
not have a control over the sum of them. For instance, it cannot be excluded 
that summing up further (or all) graphs the order of the phase transition may 
change. Thus, super daisy approximation is insufficient to solve this problem. 



4 Beyond the super daisy resummation 

To account for graphs beyond the super daisy approximation let us return to 
the Schwinger-Dyson equation in the extremum Q2"!^). Now, because in the r.h.s. 
there are contributions depending on p, the simple ansatz ( [To]) is strictly speak- 
ing insufficient. However, being interested in the infrared behavior it is anyway 
meaningful to use it at least in order to get a qualitatively correct result. After 
that we arrive at the gap equation (|23[). Here, T (|18D and T (|25|) contain the the 
higher loop graphs. In order to get an estimate of these quantities at high T we 
make the following approximations. 

Consider first T at A = 0. It is the simple 2-loop graph T(A = 0) = J ). Taking 
the contribution resulting from the zeroth Matsubara frequency only, it becomes 
T(A = 0) ~ 7T 2 with some 7 which we do not calculate explicitely. In fact, 
7 may exhibit a logarithmic dependence on T which we ignore, however. After 
that we observe that £i(g), for large q, behaves logarithmically (its superficial 
divergence degree is zero). Thereafter we approximate it by the first term of its 
hight T expansion (|57|) . It is a constant. Hence, we obtain for T ( |18D in this 
approximation 

r = 1 - n ^Tr q f3(q)^(q) ~ , (41) 

where we used ([57D and introduced the notation e = so that SAS^O) ~ e 
holds. 
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By the same procedure we obtain for f (E5J) the approximation 



Tie) 



1 - 2e 



e(2 



+ 18XTr q i3(qy-± 



\2 • 



(42) 



Here we indicated the dependence of f on e explicitely. We note Tr q j3(q) 2 = Ej(0) 
and (|5"7|) so that we obtain finally 

6e 2 (2 + e) 



r( 



1 - 2e 



:i+e) 



(43) 



This is a positive quantity of order one (note that in between T_ and T + , using 
(HD and (|0D, we have § < e < ±). 

Now we consider the gap equation in the extremum (|23|) in the leading non- 
trivial approximation and obtain 



M 



2T(e) 



m 



/ T 2 

3A — 



12 4vr 



MT m 2 M 2 



167T 2 



log 



(4vrT) s 
2m 2 



27 



+6A 2 7^ 



(44) 



This equation contains additional entries as compared to (|30|) . Due to the depen- 
dence on the mass via e it is not longer a quadratic but a higher order algebraic 
equation. Nevertheless, it is useful to rewrite it in analogy to (|31|). We obtain 



M 



3ATlYe) 



M n, 
Max 



T( 



( 3AT 



2 



da : 



7 ' 



(45) 



with 5 7 = 1 — 24A7. Now by the same arguments as above we obtain instead of 
(|34]) for the spinodal temperatures 



71 



7\ 



5i _ 2m / 



2m 
2m 5i 



3A 



167T 2 



+ 12A7 

2m ( 3A 



16tt 2 



+ 12A 7 + 



9Af(e)' 
16vr 2 



(46) 



Both these temperatures depend on 7 and T + depends in addition on f . Their 
ratio is 

T Q\ 

It remains to get some information on e. As it comes in from T + , we consider 
M Min as give by equation fl4"5D at T = T + . Taking into account that we need 
Mmiu in the leading order (~ vA) only, we obtain the equation 



M 



37Am f(e) 
2n 

10 



(at T + the square root in rhs. of Eq. (|4"5| ) is zero). In the same approximation 
we have e = and the equation can be rewritten as 



1 




(48) 



2e 



This is a fourth order algebraic equation with a solution e* = 0.318 (one of the 
other solutions is negative, two are complex). We note T(e*) = 1.57. So the gap 
( [47|) between the two spinodal temperatures becomes larger by about 50%. The 
other quantities characterizing the phase transition can be calculated in a similar 
way at last numerically. It is clear that they become all changed by a certain 
amount keeping the qualitative features unchanged. 

So we arrive at the conclusion that the higher loop graphs given by Eq. (|9|) 
in comparison to (§) do not change the order of the phase transition, making it a 
bit stronger first order. As a consequence, they do not change the problem with 
the expansion parameter discussed at the end of the preceding section. 

5 Discussion 

Let us first discuss the main features of our analysis and the results obtained. As 
it has been realized many years ago || [I4J in field theory the usual perturbation 
theory in coupling constant is not applicable in the infrared region and its series 
must be summed up to obtain more reliable results. A powerful method of re- 
summation is the second Legendre transform giving the possibility to express the 
perturbative expansion in terms of two particle irreducible diagrams. Therefore it 
accounts for the infrared divergencies coming from propagator insertions. In the 
present paper it has been applied to investigate the temperature phase transition 
in the scalar <ft 4 (x) theory with spontaneous symmetry breaking. 

As an important step, the SD equation in the minimum of the free energy 
functional was used. This considerably simplifies computations making them 
more transparent and tractable. In order to analyze the phase transition we 
adopted the ansatz (|T5| ) for the full propagator which is natural at high temper- 
ature. As we have seen in Sect. [|, all daisy and super daisy diagrams can be 
summed up by taking the 'eight '-diagram in W2, (11)- In fact, this resummation 
renders the free energy in the extremum real solving the problem of eliminating 
the instability connected with the imaginary part which constitutes a necessary 
consistency condition. The type of the phase transition was found to be of first 
order. At the same time, as it also follows from the analysis in Sec. |3], near the 
phase transition the effective expansion parameter became of order one making 
the contribution of subsequent graphs questionable. 

To go beyond the super daisy approximation we have included into considera- 
tion an infinite set of diagrams in the 2PI functional (0) and, as a consequence, the 
corresponding momentum dependent diagrams in the SD equation (P^). Their 





11 



choice is motivated by the possibility to handle them. Besides, in an extension 
to the 0(N)-model these graphs would constitute the first non-leading order for 
large N. In Sec. |] we calculated their contribution to the SD equation in the 
extremum of the free energy. Within this approximation they turned out not to 
change things essentially. Their role reduces to a redefinition of the parameters 
of the phase transition by a finite factor of order of 3/2 making the transition a 
bit stronger first order. As a consequence, the estimation of the effective expan- 
sion parameter given at the end of Sec. |3] on the super daisy resummation level, 
remains valid on the higher loop resummation of Sec. f| too. We come to the 
conclusion that even this resummation is still insufficient to get a reliable result 
from perturbation theory near the phase transition. 

We believe, anyway, that the technical tools developed here, mainly the join- 
ing of the SD equation with the necessary condition for an extremum will find 
applications in gauge theories. There, as is known, the free energy is independent 
of the gauge fixing parameter only in its minima. Therefore, additional vertex 
resummations, such as for instance done in Ref. [13 in order to obtain a gauge 
invariant result, may be avoided eventually. The realization of these ideas in 
gauge theories is a problem left for the future. 
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Appendix 



The function V\{M) = |Trln/3 appearing in fl3"9"|) reads 

K(M) = --TrlnA = T £ / i n (( 2v rZT) 2 + P + M 2 ) . 

Removing the ultraviolet divergencies the explicit expression becomes 



V X {M) 



4m 2 M 2 + M 4 



64vr 2 



In 



M 2 
2^2 



M 2 T 2 (M 
~^ b2 Vf 



(49) 



(50) 



where the last term is the temperature dependent contribution. The function Si 
can be represented as a fast converging sum over the modified Bessel function K 2 
and by an integral representation as well: 



00 1 1 

S ^ x ) = E — K ^ nx ) 



n=l 



°° , (n 2 -x 2 f 2 

, , an 

3x 2 Jx e n - 1 



(51) 
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We use 

dVl n a d2Vl n 

— = and — — = 

OV \v=m/V\ OV z \v=m/V\ 

as normalization conditions at T = 0. These conditions are chosen in a way that 
the minimum of the free energy is at T = the same as on the tree level. It is 
possible to do so as the counter-terms can be chosen at T = 0, see for example a 
recent new proof ||16|| . 

The expansion of ^(x) for small x reads 

n r \ 71-4 71-2 nx %2 f n 3 „. X \ o. 

S ^ X) = 45? ~ 12 + ~6~ + 32 ( 2 5r) + °<* > ■ (52) 

giving rise to the high temperature expansion of V\ which is at once the expansion 
for small M: 

. , -tt 2 T 4 M 2 T 2 M 3 T . . 

VAM) = + 53 

H 1 90 24 12tt v ; 



1 



64vr 2 



M 



\2 



2m 2 



Here 7 is the Euler constant. 

The function E defined in Eq. ([IB]) can be obtained by differentiating V\ 

() V 1 (M). (54) 



dM 2 

Its high temperature expansion is 

Eo(M ) = H-Ml + 'UL^ll- 2 7 ) + 2m4 + MT O . 
ov ; 12 4tt 16tt 2 \ \ 2m 2 'J J V T / 

(55) 

Similar we obtain for Ei(p) ( p0[) at p = 

E ' (0) = -aF? E » (56) 

with the high temperature expansion 

S ' (0, = ^M +O{1) - (57) 
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